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ABSTRACT 

. . . Extremely high velocity (EHV) wings, with full widths of 72 to 140 km s^^, are seen on the 

^ ! CO J ^ 3^2 lines toward W3 IRS 5, GL 490, NGC 2071, W28 A2 (G05.89-0.39), GL 2591, 

Q>^ I S140, and Cepheus A. Observations of ^^CO and ^'^CO J — 3— *-2 and J ~ 2^1 lines indicate 

that optical depth generally decreases with increasing velocity separation from the ambient 
cloud velocity. Maps of the extremely high velocity ( | — K) I ^ 20 km s~ ^ ) and the high- velocity 
(5 ^ I — Vol ^ 20 km s^^) CO emission components show that the morphology of the two 
components is similar in W3 IRS 5 and W28 A2 but may be different in GL 2591, S140, and 
Cepheus A. 

The results of our survey suggest that EHV wings are common around infrared sources of 
moderate to high luminosity [500 to (4 x 10^) Lq] in dense regions. Line ratios imply that 
the EHV gas is usually optically thin and warm. Characteristic velocities range from 20 to 40 
' km s~^, yielding timescales of 1600-4200 yr. Since most sources in this study are producing 

some ionizing photons, these short timescales suggest that neutral winds coexist with ionizing 
I photons. 

0^ ' We examined two possible sources for the extremely high velocity CO emission: a neutral 

stellar wind; and swept-up or entrained molecular gas. Neither can be ruled out. If the 
JpH| high-velocity (HV) gas is swept up by a momentum-conserving stellar wind traced by the 

Q ' extremely high velocity CO emission, most of the C in the winds from luminous objects cannot 

be in CO. If the EHV and HV forces are equal, the fraction of C in a form other than CO 
increases with source luminosity and with the production rate of ionizing photons. This trend is 
natural in the stellar wind hypothesis, but models of winds around such luminous objects are 
needed. We consider other possible chemical states for the carbon in the stellar wind. 
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1. INTRODUCTION 

Molecular outflows are common in regions of star formation (about 200 known: Lada 1985; Fukui et al. 
1993: Bachiller & Gomez-Gonzalez 1992). High-velocity (hereafter HV) molecular flows — with velocities 
ranging from a few km to about 20 km — have been studied extensively (Bally & Lada 1983; 
Plambeck, Snell, & Loren 1983; Snell et al. 1984; Richardson et al. 1985; Phillips et al. 1988). These 
molecular outflows most likely consist of ambient cloud material swept up by a faster stellar wind (Snell, 
Loren, & Plambeck 1980). 

Searches for an ionized fast stellar wind to drive the HV gas found that the momentum in the ionized 
gas was inadequate to sweep up the observed amount of molecular material (Rodriguez & Canto 1983; 
Evans et al. 1987), leading to the suggestion that the stellar wind was primarily neutral (Natta et al. 1988). 
The neutral wind gained support from the discovery of very wide wings on the spectra of H I emission 
(Lizano et al. 1988; Giovanardi et al. 1992) toward several young stars. Lizano et al. (1988), Koo (1989, 
1990), and Margulis & Snell (1989) also recognized extremely high velocity (hereafter EHV) wings on the 
CO spectra toward a few sources, suggesting that there might be CO in the stellar wind itself. In recent 
years, there have been additional observations of EHV CO outflows (Bachiller & Cernicharo 1990; Masson, 
Mundy, & Keene 1990; Bachiller, Martin-Pintado, & Planesas 1991; Mitchell & Hasegawa 1991; Chernin 
& Masson 1992; Richer, Hills, & Padman 1992). Stabler (1993) has suggested that EHV outflows are the 
manifestation in young objects of the optical jets seen in later phases of evolution. 

Calculations of the chemistry in stellar winds from low-mass stars indicated that CO could form in 
the neutral wind (Glassgold, Mamon, & Huggins 1989). Later calculations of the thermal structure of 
these winds (Ruden, Glassgold, & Shu 1990) agreed with that conclusion, and new chemical models which 
incorporated the thermal structure found that, for mass-loss rates M > 3 x 10~^ Mq yr~^ virtually all of 
the carbon is in CO ([CO]/[H] = 4 x 10"''), even though the hydrogen remains mostly atomic, unless M > 
10-'^ Mq yr-i (Glassgold, Mamon, & Huggins 1991). 

The thermal models of Ruden et al. (1990) raised questions, however, about whether a wind could 
maintain the high temperature {Tk ~ 100 K) and density (n > 10'' cm""^) deduced by Koo (1989) for 
HH 7-11, on large enough scales (-^0.1 pc) to be consistent with the observations. In addition, Masson et al. 
(1990) and Bachiller & Cernicharo (1990) found relatively narrow {AV w 20 km s"-') spectral features at 
extremely high velocities {\V — Vb| « 150 km s"'^), displaced from the center of the outflow in HH 7-11, 
suggesting that the EHV gas consists of discrete blobs, rather than a continuous wind. Similar discrete EHV 
spectral features exist in several other sources (Bachiller et al. 1990, 1991; Bachiller & Gomez- Gonzalez 
1992). The blobs may be related to Herbig-Haro objects and to localized EHV spectral features often found 
in H2O maser observations of compact H II regions (Genzel & Downes 1977, 1979). 

The nature of the EHV CO emission is currently unclear. It could arise in the stellar wind, or it could 
be ambient molecular gas which has been entrained or swept up (see Bachiller & Gomez-Gonzalez 1992 and 
Stabler 1993 for excellent reviews). 

While observing the CS J = 7— >6 line toward several sources, we noticed EHV wings on the CO 
J = 3^2 lines, which appear in the opposite sideband. We followed up by mapping the EHV gas in CO 
J = 3— >2 emission and by obtaining data on other transitions. In this paper, we report observations of 
W3 IRS 5, GL 490, NGC 2071, W28 A2, GL 2591, S140, and Cepheus A. Most previous detections of 
EHV wings have been toward sources of low to moderate luminosity {L < 1000 Lq); this work extends the 
phenomenon to higher luminosity sources. These sources, some of which are generating ionizing photons, 
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may have different wind properties if ionization plays a role in the wind generation process or in the 
chemical composition of the wind. 

We describe our observations in § 2. Section 3 explains the analysis we used to get physical parameters 
from the observations (optical depth, excitation temperature, mass, etc.). We summarize our analysis and 
discuss its implications in § 4. 



2. OBSERVATIONS AND DATA REDUCTION 

We made all observations with the 10.4 m telescope of the Caltech Submillimeter Observatory (CSO)^ 
at Mauna Kea during several runs between 1989 June and 1992 February. At the central position, all 
sources were observed in four molecular transitions: CO J = 3—^2, CO J = 2^1, ^'^CO J = 3— >2, and 
"'^^CO J = 2^1. Since most of the sources show much brighter CO J = 3^2 wings than CO J = 2— >1 
wings, they were mapped in the CO J — 3— >2 line. Table 1 summarizes the observational parameters. The 
data were obtained using the 230 GHz and 345 GHz SIS receivers in double-sideband mode with a 500 MHz 
AOS spectrometer as the backend. Telescope pointing was checked by observing CRL 618, CRL 2688, or 
Jupiter regularly and was consistent to within ^5". 

The data were taken by position switching. Table 2 lists the central and the off positions. For each 
source, an effort was made to acquire at least some of the data with an off position which was free of 
CO emission in the line core. Since some of these off positions were quite distant, other data were taken 
with more nearby off positions to improve the baselines. Therefore, the line cores in the figures may be 
affected by emission in the off positions. Some sources were observed repeatedly over several runs and the 
intensities were consistent to within ^^10%. 

The antenna temperature (T^) was calibrated by the standard chopper-wheel method, which 
automatically corrected, to first order, for the effects of atmospheric attenuation. The radiation 
temperature, Tr, following the notation of Kutner & Ulich (1981), was obtained by dividing T4 by the 
telescope main beam efficiency (rymb = Vb/Viss, where 7/h is the beam efhciency defined by Kraus (1966) 
and 77rss is the rear spillover and scattering efficiency; see Table 1). We determined rjb by observations of 
planets, primarily Jupiter (Mangum 1993). Thus, the calibration is appropriate for sources which fill the 
main beam. For sources smaller than the main beam, filling factors must be considered. Since the beam 
was about 30" for the J = 2^1 data and 20" for the J = 3^2 data, a point source centered in both beams 
will have a filling factor for the J = 3^2 data which is (3/2)^ times as high as that for the J = 2^1 data. 
For the CO J — 3^2 observations, the CS J = 7— >6 appears in the image sideband. The velocity of the 
telescope with respect to the LSR was different on each run; as a result, the intensities and shapes of the 
CS J — 7— >6 lines in Figures la-7o are not correct, and in some spectra, the CS line appears twice. 

The data were averaged for each line and a baseline was removed, using the velocity ranges listed under 
Baseline in Table 3. For the CO J = 3— >2 line, the region of the CS J = 7— >6 line was also excluded from 
the baseline fit. At the central positions, where we take line ratios, first order baselines were used. The 
resulting RMS noise levels in the baselines are given in Table 3 and spectra are shown in Figures la-7a. 



■^The CSO is operated by the California Institute of Technology under funding from the National Science Foundation, 
contract AST 90-15755. 
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Both a full-scale plot and a blowup arc shown to indicate the range and quality of the baselines. Prom the 
blown-up spectra, one can sec that the full widths to where the wing vanishes into the noise range from 
72 km s-i for W3 IRS 5 to 140 km s^^ for W28 A2 (Table 4). 

The mapping data were taken with shorter integration times, so their quality is lower. First-order 
baselines were also used for most of these spectra, but a few off-center positions required baselines up to 
third order. 

Since we are interested in the differences between HV components and EHV components of outflows, 
we have defined the inner and outer wings by examination of the CO J = 3^2 spectra. For about half 
of our sources, the boundary between the inner and outer wings was set by noticing a distinct change in 
the CO J = 3— >2 line profile, marked by the appearance of a component which declines more slowly with 
velocity away from the line core. For sources without a distinct change of slope, we set the boundary 
between our inner and outer wings at the largest velocity at which wings were previously known. For 
W28 A2, the boundaries between inner and outer wings are quite subjective because it is difficult to find a 
distinct change of slope and there is no reference useful to define the boundary. The average velocity offset 
from the line center for the boimdary between inner and outer wings was 21 km s^^, corresponding roughly 
to the boundaries chosen in other studies (e.g., Koo 1990) making it a reasonable, if not precise, choice for 
separating HV and EHV emission. The outer boundaries of the outer wings were set by our subjective 
assessment of the extremes of detectable wings. The wing boundaries are listed in Table 4 and marked in 
the blown-up spectra and the line ratio plots by dashed vertical lines. 

After defining the wing boundaries, we integrated the intensity between the boundaries and plotted 
the results as contour maps. We always set the lowest contour level higher than twice the RMS noise. 

To assess the optical depths in the wings, we examined various ratios of lines as a function of velocity. 
We resampled the spectra toward the central positions to a resolution of 2 km s~^ and formed the following 
line ratios: ^'^R32, CO J = 3^2 line divided by CO J = 2^1; ""iig, CO J = 3^2 divided by ^^CO 
J = 3^2; and "°i?2, CO J = 2->l divided by ^^CO J = 2^1. The results (Figs. Ic-lc) were blanked 
from the plot beyond the velocity where either the numerator or the denominator fell below the threshold 
(3 times the RMS noise in the resampled spectra) in Table 3. 



3. ANALYSIS 

In this section, we describe how we determined optical depths, excitation temperatures, and column 
densities from the line ratios. We also explain our calculations of the mass, momentum, and kinetic energy 
in the outflowing gas. We then compare the driving force in the EHV gas to that in the HV gas, compute 
timescales for the outflows, and collect data on the stellar luminosity and the production rate of ionizing 
photons. 



3.1. Optical Depth 
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From the isotopic line ratio, '^^"Ra, we calculated T32, the optical depth of the CO J = 3— >2 line in 
each channel, from 

isoD Trs_2 ^ l-exp(-r32) , . 

' Tfl(i3)3-2 l-exp(-r32/X)' ' 

where X (listed in Table 5) is the isotopic abundance ratio, [CO]/["'^^CO]. We assumed that X is a linear 
function of distance from the Galactic center (Langer & Penzias 1990) except for NGC 2071, where we 
assumed X = 73, the average for Ori A. The low strength in the outer wings of ^"^CO lines means that most 
sources do not have well-defined '^^°Rs at extremely high velocities. 

One way to constrain the optical depth at extremely high velocities is through the ratio of the CO 
J = 3— >2 and J = 2^1 lines; the ratio is given by 

12„ _ ^^3-2 _ [J.M.) - ■/.:i2(rbk)] [1 - exi>(-r32)] . . 

Th2-i [J.2i(rex)- J.2i(rbk)] [l-exp(-r2i)]' ^' 

where (p is the ratio of beam filling factors which accounts for differences in telescope beam sizes, Ju{T) is 
the Planck function in temperature units, T^x is the excitation temperature (assuming that the two lines 
have the same excitation temperature), and T^k = 2.7 K. In the limits of r ^ and T^x 00, ^^i?32 
approaches (i'32/^'2i)^ = (3/2)^ (for an extended source) or (^32/^21)^ = (3/2)^ (for a point source centered 
in both beams) as shown in Figure 10 (Appendix A). 

Most of the sources have '*°i?3 ratios which increase away from the line center and ^^R-i2 ratios which 
imply low optical depth in the outer wings. We assumed that the outer wings were optically thin, except 
for the red wing of S140 and both wings of GL 490, for which ^^i?32 is relatively low. For those wings, 
we assumed r32 = 2 (which is chosen arbitrarily) for all velocities. For the inner wings, we fitted T32, 
determined from **°i?3, with an exponential function (T32 = aexp[/3(F — Vq)], where a and /? are constants 
determined by a least-squares fit) for each side of the spectrum (Figs. lc-7c). 



3.2. Excitation Temperature 



For most of the outer wings, the optical depth is low and we can calculate T^x with ^^i?32 averaged 
over the wing intervals from 

12_^^2 _ ^ [Jt^32(rcx) - ^1/32(71^)] T32 , j.^^ 
[Jv2l{TeyL) — >/iy2l(Tbk)] T21 

assuming LTE, 

732 ^ 3 1 - exp(-16.597/rex) 

T21 2 exp(11.065/Tex) - 1 ' ^ ' 

It is hard to determine but we can calculate a very firm lower limit to T^x by assuming a point source 
centered in both beams [i.e., <^ = (3/2)^]. Although some maps show extended sources, the EHV wind may 
have a complex structure in position-velocity space that does not allow us to use extended emission to rule 
out a large value of 

Table 5 shows the lower limits to Te,x calculated from equation (3); they range from 16 to 94 K, 
indicating that at least some of the EHV gas is warm. If </> is closer to unity, the temperatures must be 
much higher. For example, if ^ = 3/2, an intermediate value between that expected for a point source and 
the value for an extended source, Tex ranges from 30 K to > 250 K. Any mixture of cooler gas will produce 
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a lower ratio; for example, applying this method to mixtures of two temperature components, with equal 
optical depth, would yield an average of the two temperatures. The high observed ratios are thus a strong 
indication that the bulk of the EHV gas is quite warm. 

Equation (3) cannot be used to calculate the excitation temperature in the inner wings because 
they are not, in general, optically thin. An alternative procedure would be to calculate the excitation 
temperature from the ratio of the two optical depths deduced from **°i?3 and ^-"^"R-z (for example, Snell et al. 
1984; Margulis & Lada 1985). While this method reduces systematic uncertainties due to beam filling, 
it can seriously underestimate T^^ if components with different optical depths are present (Appendix B). 
Consequently, we have not estimated T^x in the inner wings, but we show in the next section that knowledge 
of Tf-x is not essential for determining column densities. 



3.3. Column Density 



The column density of CO for each channel, Np'^ was calculated from the intensity of the CO J = 3^2 
line from 

ATf o = 1.10 X 10^5 ^zl^ Z3H (5) 

D{n,TK) l-exp(-T32) ^ ' 

where 

D{n, Tk) = /2 [J. (Tex) - J.(Tbk)] [1 - exp(-16.597/Tex)], (6) 

is the channel width in km s~^, and /2 is the fraction of CO molecules in the J = 2 state. Np*^ is in 
units of cm~^ and AV in km s"'^. The optical depth, T32, was assigned as described in § 3.1. Assuming 
that the optical depth is known, the main uncertainty is the value of the function D{n, Tk)- If all the levels 
are not in LTE, D{n,TK) depends on the density, as well as the temperature, through the factor /2. The 
LTE assumption can lead to large overestimates of N*^*^ for high temperature and modest density. 

Although Tex is not well known, LVG simulations show that D{n,TK) does not vary much within 
reasonable intervals of density, n, and kinetic temperature, Tk (Fig. 8). Therefore, we assumed D{n,TK) 
= 1.5 in our calculations, which is correct to less than a factor of 2 for 10 < Tk < 200 K and lO** <n < 10^ 
cm~^. Column densities derived from J = 3— >2 line are somewhat less uncertain than those derived from 
the J = 1— >0 line because Z)(n, Tr:) varies from 0.074 to 0.76 for the J = 1^0 transition in the same 
intervals of n and Tk- 



3.4. Mass per Channel: Two Possibilities 



The mass per channel, Mj, can be computed from the column density per channel, iV^ , the distance, 
d (Table 2), and the solid angle subtended by emission in that channel, O, by 

M,=MmHd^^^iVfoH^, (7) 

where fi = 1.3, mu is the mass of a hydrogen atom, and [H]/[C] = 2.5 x 10^ (Crevasse et al. 1991). The 
final factor, [C]/[CO], is about 8 in ambient molecular gas (Dickman 1978; van Dishoeck et al. 1992), 
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with the missing carbon presumably tied up in solid matter. For stellar winds from low-mass young stars, 
calculations indicate that molecule formation precedes dust grain formation, so that all C is in CO for 
M> 10^^ Mq yr^^ (Glassgold et al. 1991). If correct, these models say that a given N^'-' would imply 8 
times less matter if the CO emission arises from a wind than if it arises from swept-up gas. Because the 
[C]/[CO] ratio in winds from more massive stars is not well known, we will begin by assuming that carbon 
is fully associated into CO. 

In calculating mass, momentum, and energy in the next section, we give two different results for the 
EHV gas (Table 5), one for a stellar wind, with [C]/[CO] = 1, and one for swept-up gas, with [C]/[CO] = 8. 
Since the HV gas is always assumed to be swept up, we give only one result for it. 



3.5. Mass, Momentum, and Kinetic Energy 

The mass, momentum, and kinetic energy of each outflow were calculated from 

M = Y^Mi, (8) 

i 

p = Y,Mi\Vi-yo\, (9) 

i 

and 

Ek=Y,\m,{V,-Vo)\ (10) 

i 

with summations over the velocity channels, i, within each wing interval; Vo is the central velocity (see 
Table 4) of the ^^CO J = 2^1 line (W3 IRS 5 and Cepheus A) or "CO J = 3^2 line (ah other sources). 
Since we made no corrections to the velocities for projection effects, the momenta and kinetic energies in 
Table 5 are lower limits. 

Our estimates agree reasonably well with previous work. The total momentum in the wings of 

W3 IRS 5 is about 80% that found by Mitchell, Hasegawa, & Schella (1992), who included some low 
velocities in their red wing that we do not count as wing emission. Because the GL 490 maps of Snell et al. 
(1984) and Mitchell et al. (1992) cover a larger area than ours, their estimates of the mass, momentum, 
and kinetic energy of the HV outflow from CO J = 2— >1 and J = 1^0 observations are about 3-4 times 
larger than ours. The estimates of the masses, momenta, and kinetic energies of NGC 2071 by Chcrnin & 
Masson (1992) are larger (by a factor of about 5 for the outer wings and a factor of about 2 for the inner 
wings), mostly because their map is larger than ours. The W28 A2 outflow is much more energetic than the 
others in this paper. Our estimates of the momentum of the outflow agree with Harvey & Forveille (1988), 
considering the slightly different definitions of the wings and the fact that we do not include the inner red 
wing. The blue outflows of GL 2591 are much more energetic and massive than the red outflows. Mitchell 
et al. (1992) also found this asymmetry between the blue and the red outflows. The mass of the combined 
rod and blue HV outflows is about twice the mass calculated by Mitchell et al. (1992) when their values are 
corrected to our distance. Our estimates of the mass, momentum, and kinetic energy of the HV outflows of 
S140 are consistent with those of Snell et al. (1984), but exceed those of Hayashi et al. (1987), who did not 
see the full extent of even the HV wings. Our estimate of the mass in the HV flow of Cepheus A is about 
2.5 times that of Ho, Moran, & Rodriguez (1982). 
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3.6. Driving Force 

In § 4, wc will consider the possibility that the EHV wing comes from CO in the actual stellar wind. 
If the stellar wind drives the HV flow, the average driving force in the wind must equal the average force 
applied to the HV gas. The force is calculated from 

P P2 
^=7 = MS' 

where t = R/VcAi is the characteristic timescale for the outflow, R is the size of the outflow, and Vch = P/M 

is the characteristic (mass-weighted) velocity. 

The ratio of driving force in the EHV wing to the average force required to sweep up the HV wing was 
calculated from 

-Fhv Rr Rm 

where Rp, Rr, and Rm are ratios of momentum, size, and mass, respectively. It is hard to determine 
Rr from the maps for most of our sources. Since our maps and maps from previous work (see references 
in Table 6) show that the spatial extents of the EHV and the HV outflows arc similar (within a factor of 
~2), we assumed Rr 1 in our calculations. This is not true for GL 2591, which shows a very compact 
EHV source along with extended HV outflows. So the correct Rp of GL 2591 should be larger than the 
value listed in Table 6. Because the correction factor for projection should be similar for the HV and EHV 
outflows, ratios of quantities for the EHV and HV outflows should be unaffected by projection effects. 

Table 6 shows Rp for all of our sources and for some EHV outflows observed by others. Unlike our 
classification, the wings of HH 7-11 (Koo 1990) and IRAS 03282+3035 (Bachiller et al. 1991) arc divided 
into HV, IHV (intermediate high velocity), and EHV wings. For these two sources, we show two values in 
Table 6, {Fehv + Fihv)/Fhv and Fehv/{Fihv + Fhv)- Also NGC 7538 IRS 9 (Mitchefl & Hasegawa 
1991) has an EHV wing detected only to the red side of the spectrum. Again, we list two values in Table 6, 

FpedEHV / FpedHV and FpedEHV / {FRedHV + FbIucHv)- 

Table 6 also shows the characteristic timescale for the EHV outflows t = R/Vch- When the map shows 
a clear separation between blue and red lobes, the size of the outflow, R, is calculated from the distance 
between red and blue peaks. If the blue lobe is superposed on the red lobe, we take the size of the FWHM 
contour as R. Projection effects, not corrected for in our calculation, can increase (if flow is in the plane of 
the sky) or decrease (if flow is along the line of sight) the timescale by substantial factors. In addition, t 
would be smaller if we used the maximum velocity, as is often done, rather than Vch- On the other hand, 
Parker, Padman, & Scott (1991) have argued on the basis of detection statistics that timescales derived 
from this method must underestimate the true ages of outflows, by an order of magnitude on average. 



3.7. Other Information 

Since the discussion in § 4 will need the luminosity and production rate of ionizing photons, these are 
given in Table 6. The luminosity was taken from the literature, corrected if necessary to the distances we 
have adopted. The rate of production of Lyman-continuum photons, Nl was calculated from 

^43"0-irf25^T-0-45^ (13) 
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where Ni, is in units of s""'^, v the frequency in GHz, d the distance in kpc, 5*^ the radio continuum flux 
in mJy, and the electron temperature in units of 10** K (Rubin 1968). We assumed Te = 8000 K and 
used the radio data given in Table 6. To avoid optical depth problems, we used the highest frequency radio 
continuum data available which was not judged to be contaminated by dust emission. 

In some of these sources, the radio continuum emission may arise in an ionized stellar wind. In these 
cases, Nl calculated using equation (13) will overestimate the production of Lyman-continuum photons 
because Balmer-continuum photons may also ionize the dense winds (Thompson 1984). We have computed 
the mass-loss rates in ionized winds {Mionized) by assuming that all the radio emission in Table 6 comes 
from an ionized wind. We have used equation (30) of Panagia (1991), appropriate for a spherical, isothermal 
wind, and assumed a terminal wind velocity of 300 km s^^ and an electron temperature of 10^ K. If the 
wind is confined to a narrow jet, this equation will overestimate the mass-loss rate (Reynolds 1986). In 
addition, any contribution to S^, from an optically thin H II region will cause us to overestimate Mionized- 

For comparison, we have also calculated the mass-loss rate in the stellar wind based on the HV CO 
emission, assuming momentum conservation, using 

Mco = 7^, (14) 

tVw 

where P and t are the momentum and timescale for the HV flow, and Vyj is the stellar wind velocity, 

assumed to be 300 km s"-'^. Mionized and A^co are given in Table 6. In cases where our maps are incomplete, 
Mco will be underestimated, but if Parker et al. (1991) are right about the ages of outflows, Mco based 
on the usual method will be overestimated. In defense of the traditional method, we note that Natta & 
Giovannardi (1990) derived mass loss rates from Na I lines for several stars which agreed very well with 
rates derived from CO. 



4. CONCLUSIONS 



In this section, we will summarize the results and then consider diflferent ideas for the origin of the 
EHV gas and the nature of the stellar wind. 



4.1. Summary of Results 



Our detection of EHV wings was a result of serendipity while conducting a CS study. The selection 
criteria involved previous detections of dense gas rather than any knowledge of molecular outflows, although 
all the sources were known to have HV outflows. Thus, it is probably signiflcant that almost all the sources 
in the study did show EHV wings. EHV outflows are apparently common in regions similar to those in our 
sample, dense regions of moderate to high mass star formation, with luminosities ranging from 500 Lq to 

4 X 10^ Lq. 

The J = 3^2 line of CO shows the EHV wings more clearly than lower J lines, indicating that it is a 
good way to search for EHV gas, given an appropriate receiver and site. The high ratios of the J = 3— >2 
to J = 2^1 EHV wing emission indicate that the EHV gas is generally (GL 490 is an exception) optically 
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thin and warm {Tex > 20-90 K). While the formal limits on Tex in Table 5 are not very high, these limits 
assume the highest possible value for the ratio of filling factors. Lower values of this ratio would produce 
much higher Tex, comparable to the temperatures needed to excite the CO J = 14^13 and J = 16^15 
lines observed toward many of these sources (Stacey 1992), and impossibly high in some cases. Howe et al. 
(1987) have detected J = 14^13 line of CO toward NGC 2071. They remark that the warm component 
responsible for the far-infrared CO emission can account for the high- velocity gas seen in the millimeter CO 
lines if the density is near 5 x 10^ cm~^ and the temperature of the gas lies between 250 and 750 K. 

The spectra indicate molecular gas moving at characteristic velocities (in the sense of momentum over 
mass) of at least 20-40 km s~^ with respect to the cloud rest velocity, without projection corrections. 
It is hard to accelerate molecular gas to such high velocities without dissociation. Based on the size 
of the outflows, these velocities result in characteristic times of 1600-4200 yr. Thus, these stars have 
been producing stellar winds, probably neutral winds, within quite recent times, even though many now 
are producing enough ionizing photons to have substantial H II regions. One explanation, particularly 
attractive in distant, very obscured regions, is that the source driving the outflow is separate from the 
source producing the ionizing photons. In the case of S140, however, the evidence points to IRS 1 as the 
source of the neutral wind, even though it may be a main sequence star (Evans et al. 1989). 

In contrast to recent spectroscopy of nearby, less luminous sources with EHV gas, our spectra are 
dominated by fairly smooth profiles, with discrete velocity features being rare; in addition the maps show 

clumpy emission only in the case of S140. On the other hand, our maps are often incomplete or not 
completely sampled. The weakness of the EHV wings, together with their likely high Tg^, suggests a low 
filling factor, which could result from many small clumps. In this case, the smooth profiles could just be a 
result of inadequate spatial resolution. Another shared spectral characteristic is the presence of distinct 
changes of slope between the EHV and HV wings, with the slope discontinuities in NGC 2071 and GL 2591 
being the most dramatic. 



4.2. Nature of the EHV Gas 

We consider a quantitative test of the hypothesis that the EHV gas traces the stellar wind. Giovanardi 
et al. (1992) have shown that the H I winds in L1551 and HH 7-11 are sufficient to drive the HV molecular 
outflow by momentum conservation. We apply a similar test to the EHV gas, traced in our case by 
CO. In Table 6, we give Rp, the ratio of the force available in the EHV gas to that required to produce 
the momentum in the HV gas, averaged over the lifetime of the HV flow. We also include Rp for some 
previously studied sources. 

While the ratios for previously studied sources are consistent with unity within likely uncertainties, the 
ratios for our sources are not. The largest ratio found in our observations is 0.10 (GL 490) and the average 
is 0.045. Low ratios appear to rule out the stellar wind interpretation for our sources, but there are several 
reasons to suspend judgement. Since the EHV gas is spread over a large range of velocities, sensitivity 
limitations could cause an underestimate of the emission. Second, our division between HV and EHV gas 
is not precise in sources without clear spectral breaks; moving the boundary between HV and EHV gas 
to lower velocities would raise the ratios. Examples of this eff'ect can be seen in the different ratios for 
HH 7-11 and IRAS 03282-f3035, depending on whether the intermediate-velocity wing is assigned to the 
HV or EHV wings. Third, by drawing our boundaries in velocity rather than fitting a series of Gaussians, 
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we are neglecting gas at low projected velocities, a more serious omission for the EHV gas than for the HV 
gas. Fourth, incomplete maps in some sources introduce uncertainties; for example, the data of Chernin 
& Masson (1992) yield Rp about 3 times our value for NGC 2071. These four factors typically introduce 
factor of 2-3 uncertainties, so are unlikely to change the conclusions. 

Perhaps more important is our assumption that all the C is associated into CO in calculating Rp for 
the stellar wind hypothesis. If, on average, only 4.5% of the C is in CO, the average Rp would be unity. 
The chemical calculations of Glassgold ct al. (1991) find CO abundances as low as [CO]/[C] — 0.04 only 
for low mass- loss rates (M ^ 3 x 10"'' M© yr"^), but our EHV outflow sources have much larger mass-loss 
rates (M > 10~* Mq yr~^). However, the calculations of Glassgold et al. (1991) were for cool (T* = 
5000 K) stars, while the stars in this study arc much more luminous and most arc producing some ionizing 
photons. Ionizing and dissociating photons from these stars may cause lower CO abundances in the wind. 
Glassgold et al. (1991) did consider a model with M = 3 x 10"^ Mq yr'^ which had aT^ = 6000 K. Even 
this modest increase in r,^ caused the CO abundance to drop by three orders of magnitude. Glassgold et al. 
(1991) also note that they have ignored dust formation in the wind which could shield molecules. 

We need models of winds from stars with L > 500 Lq and higher temperatures to determine CO 
abundances in the winds. Empirically, we can assess whether the nature of the stars may be affecting the 
CO abundance by plotting Rp versus luminosity and versus A''^, the rate of production of ionizing photons 
(Fig. 9). The plots indicate that Rp declines sharply with increasing luminosity and increasing Nl, with 
the latter trend being particularly striking. If the CO abundance decreases with ionizing photon rate, the 
EHV CO emission could still be tracing the stellar wind. W28 A2 is unusual in having an impressive EHV 
wing while producing a large N^. As a distant and complex region, W28 A2 may be a case in which the 
outflow is driven by a separate source. However, the outflow must still avoid ionization by the source of the 
ionizing photons for the extended H II region. 

We conclude that an origin for the EHV wings in molecules formed in a neutral wind cannot be clearly 
ruled out. Even if the EHV CO does arise in the wind, it is not a useful probe of the momentum in high 
luminosity sources, since the CO is unlikely to be fully associated. 

An alternative origin for the EHV gas is ambient gas which has been entrained or swept-up by the 
wind. In cases where there is a break in the slope in the profile, the process must be distinct from that 

which produces the HV wing. Chernin & Masson (1992) argue that it is difficult to accelerate molecules to 
the extremely high velocities without dissociation, so they favor a picture in which the EHV molecules have 
reformed after the passage of a fast shock. Models of such shocks indicate that the molecules reform when 
the column density of the shocked material reaches about 1-3 x 10^° cm^^ (Neufeld & Dalgarno 1989a; 
Hollenbach & McKec 1989). By = 3 x 10^" cm~^, the carbon is fully associated, with an abundance 
of 4 X lO^'', implying a CO column density of about 1 x 10^^ cm~^, and the gas has a temperature of 
about 500 K. The high temperature is consistent with observations of high- J CO emission in several of our 
sources (Howe ct al. 1987; Stacey et al. 1987) and with our high i2i?32, but our column density estimates 
for the EHV CO at the central positions range from 2 x 10^"'' to 1.2 x 10^^ cm^^, with only W28 A2 having 
the column density expected in this picture. Larger column densities may be present if the emission comes 
from many small regions, as seems likely on other grounds. The recent detection of SO in NGC 2071 clearly 
indicates that shock chemistry is active, but the spatial and velocity structure of the SO emission is hard to 
interpret simply in terms of shock models (Chernin & Masson 1993). 

In the picture of outflows proposed by Stabler (1993), the stellar wind sweeps up and entrains gas in a 
turbulent mixing process, which may allow acceleration to EHV velocities without dissociation. In this case. 
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the original molecular cloud abundances may be preserved, with only 12% of the C in CO. We can then 
compare the mass, momenta, and energy of the EHV gas to that of the HV gas, using the second columns 
for each in Table 5. It is clear that the HV gas carries most of the mass and momentum in all cases; in a 
few cases, the energy in the EHV wings is comparable to that in the HV wings, but it is usually not. Thus, 
the EHV wings, if swcpt-up matter, do not change the overall energetics of the outflows. The Rp values are 
still less than 1 (Table 6), but the ratios are now sufHciently close (0.36 in the average) that uncertainties 
could explain the difference. In the picture of Stabler, in which momentum is transferred through turbulent 
mixing, one might expect a similar force to have been applied to gas at all velocities. 

A different model of outflows, with collimated jets as the driving source, has been suggested recently 
(Masson & Chcrnin 1993; Raga ct al. 1993). .Jets arc plausible driving sources only if Parker et al. (1991) 
are correct in arguing that traditional analyses have underestimated lifetimes and overestimated Mco 
by a factor of 10. In the jet models, the EHV CO emission is produced by the bow shock, where the 
jet encounters the ambient cloud. This simple model would predict a single EHV feature at the end of 
each lobe, rather than an extended wing, but Masson & Chcrnin (1993) proposed that the jets wander, 
producing multiple EHV blobs. Whether this model can reproduce the rather smooth EHV wings and the 
large velocity extents that we observe remains to be seen. 

In general, the picture of the EHV wings as gas originally in the ambient cloud is certainly plausible. 
However, this picture provides no obvious explanation for the rather striking trend of Rp with Nl in 
Figure 9, nor for the fact that i?ir ~ 1 in some sources with low luminosity and N^. 



4.3. Where is the Ccirbon in the Wind? 

If the EHV material is, in fact, ambient material swept up by a faster stellar wind, we would still need 
to find the underlying stellar wind driving both the HV and EHV components. In young low-luminosity 
stars, ionized winds, if they are present at all, are insufficient to sweep up the molecular gas (Evans et al. 
1987). The situation for more luminous objects is less clear. In Table 6, we compare mass-loss rates inferred 
from the radio emission {Mi^nized) with the stellar wind M inferred from the HV CO (Mco) assuming 
momentum conservation (eq. [14]). Only GL 490 and radio component 2 in Cepheus A actually have the 
properties expected for spherical winds, and Mionized for those sources is an order of magnitude less than 
Mco- The estimates of Mionized are upper limits for the other sources, and Mco estimates are often lower 
limits. Still, only W28 A2 and GL 2591 have Mionized comparable to Mco- The radio source in GL 2591 
does not coincide with the infrared source, making its status unclear. W28 A2 has been classified as an 
ultracompact H II region with a shell structure (Wood & Churchwell 1989) and its radio emission is too 
strong to arise in an ionized stellar wind. 

The best interpretation of the radio continuum results for our sources is that the true value of Mionized 
is much less than Mco in all of the sources in the sample. This leads to the conclusion that the winds 
are primarily neutral, even in these more luminous objects. Detection of H I wings toward DR 21 (Russell 
et al. 1992) supports the idea that luminous stars may also have neutral winds. How a star producing a 
substantial number of ionizing photons also drives a neutral stellar wind is a puzzle worthy of theoretical 
consideration. We will confine ourselves to asking where is the carbon in the stellar wind? Since these 
are oxygen-rich objects, the carbon is unlikely to be found in solid particles, leaving CO, C°, and C+ as 
possibilities. 
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In low luminosity sources, the carbon may well be in CO and already detected since Rp ~ 1- 
Observations of H I in some low-luminosity sources over large velocity intervals indicates that there is atomic 
or molecular material at velocities up to 150 km from line center (Lizano et al. 1988; Giovanardi et al. 
1992). Could we have detected wings of this width if the gas contains enough CO to make Rp = 1? A wind 
with a rectangular profile and a terminal velocity of 150 km would produce a brightness temperature 
less than 0.2 K in the J = 3— »2 line of CO for all sources except W28 A2, where it would be more like 2 K. 
Given that ±150 km would extend to the ends of our baseline and given the systematic uncertainties 
in the baselines, it would be difficult for us to rule out the existence of this kind of CO component, even if 
the spectral shape had some degree of central peaking. Continuing improvements in instrumentation and 
observing techniques may eventually permit a meaningful search for even higher velocity molecular gas. 

Another possibility is that the C is neutral and atomic. We have estimated the strength of the EHV 
wings in '^Pi^^Pq for the sources in our sample, assuming that the gas is hotter than 50 K, all of the 
carbon not in CO is in C", and the force in the EHV wind matches that in the molecular HV wind. With 

these assumptions, it should be possible to detect C° emission from the EHV gas in at least half the sources 
in our sample. For most, the [C I] line would be detectable at some EHV velocities even if [C°]/[C] is as 
low as 0.1. Walker et al. (1993) have observed HV wings in the ^Pi^^Pq transition of C I at 492 GHz in 
several sources, including three of the sources discussed in this paper. They found that the C° abundance 
was comparable to the abundance in quiescent cloud material. Walker et al. (1993) would not have been 
able to detect EHV wings, but longer integration time under good conditions should be able to determine if 
a substantial fraction of the carbon in the stellar wind is C°. 

The carbon could also be in the form of C+, if it is ionized by stellar photons. High- velocity shocks 
could also produce substantial ultraviolet emission, but shock models for dense gas predict 0+ column 
densities which are always less than those of CO (Neufeld & Dalgarno 1989b). The strong trend for Rp 
to decrease as Nl increases may indicate a chemical shift from CO to C+, due to the stellar ultraviolet 
emission. Estimates based on assumptions similar to those used for [C I] predict strong [C II] emission from 
a mostly 0+ wind. The sensitivity of current airborne spectrometers to the 158 /itm ^^3/2— >^-Pi/2 line is 
sufficient, in principle, to place meaningful limits on the amount of 0+ in the EHV winds in our sources, 
but practical limitations would make it very hard to detect very weak, broad wings on a narrow bright line 
in the far-infrared. 

While the observations are difficult, failure to detect C in any form at the level required by the usual 
analysis of the CO HV flows would support the idea that outflow ages have been underestimated and lend 
credence to models of jet-driven outflows. 
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helpful comments. This work was supported by NSF grant AST-9017710, by a grant from the W. M. Keck 
Fomidation, by a grant from the Texas Advanced Research Program to the University of Texas, and by a 
David and Lucile Packard Foundation Fellowship. 
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To understand the general behavior of ^^i?32, we can simphfy equation (2) by assuming that the optical 
depths follow the LTE relation (eq. [4]). The largest possible value of the ratio of beam filling factors, (f> 
is (3/2)^, assuming a point source centered in both beams. The resulting ^^i?32 is plotted in Figure 10 as 
a function of the optical depth for various values of excitation temperature. For an optically thick cloud, 
"^^i?32 is quite insensitive to both t and T^x- For an optically thin cloud, '^^i?32 can be as low as 0.44 (for 
Tex = 10 K and </> = 1) and as high as (3/2)'* [for Tex oo and </> = (3/2)^] depending on the temperature 
and source extent. Therefore, ^^i?32 larger than ~(3/2)^ implies that the cloud is optically thin, but *^i?32 
smaller than ~(3/2)^ does not give any clue to the physical conditions unless the source extent is known. 

One might regard ^^i?32 as a good temperature indicator for optically thin clouds based on Figure 10, 

but it is hard to dc'tc^rminc (f>, the ratio of beam filling factors. The lower limits to T^x in Tabic 5 were 
obtained by assuming a point source. In addition, mixtures of gas at different temperatures will result in an 
estimate of T^x less than that of the hotter gas. If both components have equal optical depth, the derived 
temperature is close to the average of the two components, according to simulations that we have done. 
Simulations with two equal temperature components of differing optical depth, but both still optically thin, 
show that the true temperature will be somewhat underestimated by this method. Thus, the excitation 
temperatures we have derived from this method should be regarded as quite strong lower limits. 



B. Temperature Estimates from the CO Optical Depth Ratio 



It is in principle possible to determine the excitation temperature from the ratio of the optical depths 
of two adjacent transitions assuming LTE and, in the process, avoid some of the problems with beam filling 
factors alluded to in Appendix A. Consider two transitions, CO J+1 J (labeled by u) and CO J J-1 
(labeled by Z); then 

Tu ^ + 1 1 - cxp j-hl^u/kTcx) ^ 
n J cxp(hiyi/kTcx) - I 

Snell et al. (1984) and Margulis & Lada (1985), using this method for several outflows with CO J = 2— >1 
and CO J = 1^0 data (calculating r2i and no from '''°i?2 and using eq. [1]), found oddly low 

excitation temperatures {Tex = 6.5-25 K) except in Ori A. Although we have determined optical depths of 
CO J = 3^2 and CO J = 2— >1 lines for most of the inner wings and some of the outer wings, we did not 
use this method to calculate the excitation temperatures. The reason is that outflows may be clumpy and 
there may be several unresolved components in our main beam. 

Assuming that there are two components, A and B, in our main beam, one with t„a and the other with 
r„B, but both with the same excitation temperature, Tex, then tia and tje can be deduced by equation (Bl). 
The observed radiation temperature is roughly the sum of those from each component. The procedure 
would derive the optical depths, Tu,obs and Ti^obs, assuming all the fllling factors are about the same, from 

l-exp(-Tobs) l-exp(-rA) + l-exp(-TB) 

- (B2) 



1 - exp(-Tobs/-'«^) l-exp(-rA/X) + 1 - exp(-rB/X) ' 

where t can be either t„ or t/. 

The resulting Tex,obs (excitation temperature calculated by eq. [Bl]) is shown as a function of Tex 
(real excitation temperature of the cloud) in Figure 11, assuming two components with equal filling factor 
with one component having r = 1 and the other r = 10. These curves clearly show that this method can 
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fail badly as a temperature estimator if there are two or more components of different optical depth in the 
beam. The qualitative result is true for any pair of opacities as long as they differ by a factor of 2-3 or more 
and the component with larger opacity has r > 1. If the low temperatures derived by this method from 
the J = 2— >1 and J = 1— >0 data are an artifact of using the technique in the presence of inhomogeneities 
in the sources, Figure 11 shows that a similar analysis using the J = 3^2 and J = 2^1 transitions should 
yield higher temperatures. Indeed, the typical T^x^s we derive by this method are twice those derived in the 
earlier work with the J = 2^1 and J = 1— >0 lines. 



- 16 - 



TABLE 1 



Observational Parameters 



Telescope Transition 




V 






'7m6 






(GHz) 


(km s^"*^) 






CSO i^cO J = 2- 


>1 


220.3987 


1.37 


32" 


0.71 


CO J = 2- 


>1 


230.5380 


1.28 


30 


0.71 


i^CO J = 3- 


>2 


330.5880 


0.91 


21 


0.55 


CO J = 3- 


>2 


345.7960 


0.86 


20 


0.55 


" Velocity resolution 


in 


km s~^. 









Beam size (FWHM) in arcseconds. 
Main beam efSciency (Mangum 1993). 



TABLE 2 



Source Parameters 



Source 


ai950 


^1950 


. 

(kpc) 


References 






Line'^ 




W3 IRS 5 


02''21"53!1 


61°52'22" 


2.3 


1 


2400" 


0" 


3 2, 133-2, 132-1 














1200 





2 1 




GL 490 


03 23 39.2 


58 36 36 


1.0 


2 


-1800 





3-2, 2-1, 133-2, i32_ 


-1 


NGC 2071 


05 44 30.6 


00 20 42 


0.39 


3 


-2700 


900 


3-2, 133-2, 132-1 














-2700 





2-1 




W28 A2 


17 57 26.8 


-24 03 54 


3.0 


4 


-600 





3-2, 2-1, 133-2, 132- 


-1 


GL 2591 


20 27 35.8 


40 01 14 


2.0 


5 





10800 


3-2, 2-1, 133- 2, 132- 


-1 


S140 


22 17 41.2 


63 03 45 


0.91 


6 


-600 


-600 


3-2, 2-1, 133-2, 132- 


-1 


Cepheus A 


22 54 19.0 


61 45 47 


0.725 


7 


-780 





3-2, 2-1, 133-2 

















1200 


132-1 





" Distance from the Sun in kpc. 

* Coordinates of off positions relative to the central position. 

" 3-2: CO J = 3^2; 133-2: i3C0 J = 3^2; 2-1: CO J = 2-^1; 132-1: i3C0 J = 2-»l. 

References.— (1) Georgelin & Georgelin 1976; (2) Snell et al. 1984; (3) Anthony-Twarog 1982; 
(4) Harvey & Forveille 1988; (5) Merrill & Soifer 1974; (6) Crampton & Fisher 1974; (7) Blaauw, 
Hiltner, & Johnson 1959. 
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TABLE 3 



Line Parameters 



Source 


T • (7 

Line 




Baseline'' 


Noise (K) 


Threshold (Kj 


W3 IRS 5 


3 2 


( 150, 


90) 


(10, 40) 


0.115 


0.264 




9-1 




-90) 


('M) 1 {){)) 

1 Ovj ^ -L yjyj 1 


u.uoo 


u.uyu 




133-2 


(-150, 


-7(1"! 


I L/ 5 -L \J\J 1 


0.116 


0.285 




132-1 


(-150, 






0.019 


0.048 


GL 490 


3-2 


(-100, 


-75) 


(50, 100) 


0.056 


0.126 




9-1 


(-^ nn 


ou 1 




U.IOO 


n 9Q7 




133-2 


(-100, 


-50) 


('40 TOO"! 


0.291 


0.645 




132-1 


(-100, 


-50) 


(40 100) 


0.114 


0.264 


NGC 2071 


3-2 


(-130, 


-60) 


(100, 150) 


0.091 


0.195 




9 1 




\J\J J 


Cqn 1 i^nA 


u.uoo 


U.Ucj't 




133-2 


(-130, 


-60) 


(QO 1 501 


0.064 


0.135 




132-1 


(-130, 


-50) 


(60 1 50") 


0.035 


0.084 


W28 A2 


3-2 


(-100, 


-90) 


(180, 220) 


0.225 


0.579 




9—1 


('— 9nn 


-100) 


noo 


u.uyu 


n 998 




133-2 


(-180, 


-80) 


(70 ISO"! 


0.222 


0.546 




132-1 


(-150, 


-100) 




0.211 


0.516 


GL 2591 


3-2 


(-120, 


-50) 


(25, 85) 


0.211 


0.528 




2-1 


(-120, 


-50) 


(25, 85) 


0.070 


0.180 




133 2 


(-120, 


-50) 


(25, 85) 


0.189 


0.432 




132 1 


(-120, 


-40) 


(20, 85) 


0.026 


0.072 


S140 


3-2 


(-90, - 


60) 


(30, 80) 


0.084 


0.219 




2-1 


(-100, 


-60) 


(30, 100) 


0.046 


0.117 




133-2 


(-100, 


-60) 


(30, 100) 


0.078 


0.180 




132-1 


(-100, 


-60) 


(30, 100) 


0.037 


0.102 


Cepheus A 


3 2 


( 200, 


150) 


(80, 120) 


0.422 


0.900 




2-1 


(-150, 


100) 


(100, 150) 


0.285 


0.639 




133-2 


(-100, 


-50) 


(30, 100) 


0.120 


0.306 




132-1 


(-100, 


-50) 


(30, 100) 


0.027 


0.069 


" See Table 2. 



° Intervals considered in the baseline determination in km s~ 
RMS of noise with the resolution in Table 1. 
Thresholds for the line ratios in data with Av = 2 km s~i. 
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TABLE 4 



Wing Boundaries 



Source 




Blxje 








Red 




Full Width 




Outer 


Inner 






Inner 


Outer 


OF Wing 


W3 IRS 5 


-76 


-60" 


-50 


-37 


-28 


-18" 


-4 


72 


GL 490^^ 


-76 


-40^ 


-16 


-13 


-8 


14& 


36 


112 


NGC 2071* 


-50 


-10" 





11 


18 


28" 


68 


118 


W28 A2f 


-70 


-20^^ 


2 


9 


18 


32^= 


70 


140 


GL 2591S 


-49 


-21" 


-11 


-6 


3 


15" 


25 


74 


S140'' 


-56 


-30^ 


-12 


-7 


-4 


6" 


30 


86 


Cepheus A' 


-81 


-29" 


-17 


-12 


-5 


7" 


51 


132 



Note. — All velocities are in km s~^. 

" The boundaries between inner and outer wings were set by noticing distinct 
changes of the slope in the CO J — 3^2 spectra. 

The CO J = 3^2 spectra do not show clear change of slope, and the inner 
wings were defined from previous observations of HV outflows. 

The CO J = 3^2 spectra do not show clear change of slope, and the 
boundaries between inner and outer wings were set subjectively. 

Our inner wings coincide with the wings seen by Snell et al. (1984). The 
full widths of the EHV CO J = 1^0 Une of Koo (1989) and MarguUs & Snell 
(1989) are wider than ours by 40 and 13 km s~^, respectively, probably because of 
higher signal-to-noise ratios. 

Our inner wings coincide with the entire wings seen by Bally (1982), while 
our outer wings are wider by ~40 km s~^ on each side, similar to those of Chernin 
& Masson (1992). 

The full width of the wings is comparable to the width in CO J = 1— >0 line 
(Harvey & Forveille 1988). 

3 Our inner wings have the same boundaries as defined in Lada et al. (1984) 
and are similar to those of Mitchell, Maillard, & Hasegawa (1991). 

Our inner blue wing has the same extent as the blue wing of Hayashi et al. 
(1987) and Snell et al. (1984). The full width of the EHV CO J = 1-^0 Hne of 
Koo (1989) is wider than ours by 30 km s-\ but that of Margulis & Snell (1989) 
is narrower than ours by 11 km s~^. 

' Our inner wings correspond to the wings shown in Figs. 36-35 of Ho et al. 
(1982). Our outer blue wing is wider than the wing in Fig. 3/i of Ho et al. (1982), 
and our outer red wing was not seen by them. 
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TABLE 5 



Physical Parameters of Outflows 









ex 


IVl 


IVl 


pc 


pd 




Ei 


Source 


Wing 


A 








{Mq km 


(10^5 erg) 


Wo irvo O 


Outer blue 


(id 
DO 




U.Ui4 


nil 

U.ii 


u.oo 


1 

0. 1 


nil 

U. 11 


n R7 

U.O 1 




Inner blue 








1.5 




94 




A 




Inner red 








6.2 




78. 




10. 




Outer red 






0.026 


0.21 


0.63 


5.0 


0.15 


1.2 


GL 490^= 


Outer blue 


61 




0.0084 


0.067 






n 1 

u. ±o 


1 




Inner blue 








2.3 




1 Q 




9 1 

Z . 1 




Inner red 








1.9 




18. 




2.1 




Outer red 






0.018 


0.14 


0.60 


4.8 


0.21 


1.6 


NGC 2071"= 


Outer blue 


73 


>17 


0.0049 


0.039 




1 1 








Inner blue 








0.28 




O.-J 




yf.OO 




Inner red 








0.28 




2.7 




0.28 




Outer red 






0.0029 


0.023 


0.083 


0.67 


0.027 


0.21 


W28 A2 


Outer blue 


40 


>94 


0.16 


1.3 


7 n 




9 
o.^ 






Inner blue 








38. 




•inn 




7^ 
1 o. 




Outer red 






0.14 


1.1 




O 1 . 


1 1^ 
i.O 


1 9 
IZ. 


KjrLi ZOyi 


Outer blue 






U.UoU 


U.Z4 


0.66 


5.3 


0.15 


1.2 




Inner blue 








31. 




240. 




19. 




Inner red 
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" The isotopic abundance ratio, [CO]/[^^CO] (Langer & Penzias 1990). 
^ Lower limits of T^x calculated from eq. (3) assuming (p = (3/2)^. 
^ Assuming that the sources are neutral atomic wind. 

Assuming that the sources are swept-up molecular gas. 
^ Maps cover only the central portions of the outflows. 



TABLE 6 



Ratios of Driving Forces and Related Quantaties 
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(10 Mq yr ) 
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W3 IRS 5 
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0.15 


1 


3100 
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2 


15.4 


7.0 


3 


45.6 


1.1 


GL 490 


0.10 


0.82 


1 


2300 


1.3 


1700 


4 


23.1 


3.2 


5 


44.6 


0.14 


NGC 2071 


0.082 


0.66 


1 
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0.3 
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6 
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9.2= 


7 


44.2 


0.09 
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0.36^ 
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8 


15.0 
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9 


48.8 
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GL 2591 


0.0085 


0.068 


1 
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6.8 
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10 


5.0 
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11 


46.5 


9.1 


S140 


0.015 


0.12 
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12 


15.0 


7.7'' 


13 


44.8 


0.29 


Cepheus A 


0.041 
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1 
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14 
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5 


44.9 


0.26 
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B. Previous Observations 


HH 7-11 


1.2^ 


9.8^ 


15 






115' 


16 


5.0 


0.8'" 


7 


43.0 


0.02 


HH 7-11 


0.65'= 


5.2*= 


15 




















HH 7-11 


0.33 


2.7 


17 




















IRAS 03282+3035 


0.72^ 


5.7^ 


18 






<1.5 


18 












IRAS 03282+3035 


0.23*= 


1.8'= 


18 




















NGC 2071 


0.24 


1.9 


19 




















NGC 7538 IRS 9 
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" Ratio of driving forces assuming the EHV outflow is the stellar wind ([C]/[CO] — 1). 

^ Ratio of driving forces assuming the EHV outflow is swept-up molecular gas ([C]/[CO] — 8). 

^ Characteristic timescale of the EHV outflow. ^ Rate of production of Lyman-continuum photons. 

" Prom NGC 2071 IRS 1 alone. f From blue outflows only. 

» Prom GL 2591 (1) alone. Prom S140 IRS 1 alone. 

' Prom Cepheus A (2) alone. 

Assumes IHV outflow is part of EHV outflow. ^ Assumes IHV outflow is part of HV outflow. 

' Prom SVS 13 alone. Prom NGC 1333 (5) alone. 

" Prom red EHV and red HV outflows only. ° From blue HV, red EHV, and red HV. 

References. — (1) This paper; (2) Werner et al. 1980; (3) Colley 1980; (4) Harvey et al. 1979 (corrected according to our distance); (5) Simon et al. 
1983; (6) Butner et al. 1990; (7) Snell & Bally 1986; (8) Emerson, .Tennings, & Moorwood 1973: (9) Wood & Churchwell 1989; (10) Lada et al. 1984: 
(11) Campbell 1984; (12) Lester et al. 1986; (13) Evans et al. 1989; (14) Koppenaal et al. 1979 and Evans et al. 1981; (15) Koo 1990; (16) Cohen & 
Schwartz 1987; (17) Bachiller & Cernicharo 1990; (18) Bachiller et al. 1991; (19) Chernin & Masson 1992; (20) Mitchell & Hasegawa 1991. 



REFERENCES 



Anthony-Twarog, B. J. 1982, AJ, 871, 1213 
Bachiller, R., & Cernicharo, J. 1990, A&A, 239, 276 

Bachiller, R., Cernicharo, J., Martm-Pintado, J., Tafalla, M., & Lazareff, B. 1990, A&A, 231, 174 

Bachiller, R., & Gomez-Gonzalez, J. 1992, Astron. Astrophys. Rev., 3, 257 

Bachiller, R., Martfn-Pintado, J., & Planesas, P. 1991, A&A, 251, 639 

Bally, J. 1982, ApJ, 261, 558 

Bally, J., & Lada, C. J. 1983, ApJ, 265, 824 

Beichman, C. A., Becklin, E. E., & Wynn- Williams, C. G. 1979, ApJ, 232, L47 
Blaauw, A., Hiltner, W. A., & Johnson, H. L. 1959, ApJ, 130, 69 

Butner, H. M., Evans, N. J., II, Mundy, L. G., Natta, A., & Randich, M. S. 1990, ApJ, 364, 164 

Campbell, B. 1984, ApJ, 287, 334 

Chernin, L. M., & Masson, C. R. 1992, ApJ, 396, L35 

1993, ApJ, 403, L21 

Claussen, M. J., et al. 1984, ApJ, 285, L79 
Cohen, M., & Schwartz, R. D. 1987, ApJ, 316, 311 
CoUey, D. 1980, MNRAS, 193, 495 

Crampton, D., & Fisher, W. A. 1974, Pub. Dom. Ap. Obs., 14, 12 
Dickman, R. L. 1978, ApJS, 37, 407 

Emerson, J. P., Jennings, R. E., & Moorwood, A. F. M. 1973, ApJ, 184, 401 
Evans, N. J., 11, et al. 1981, ApJ, 244, 115 

Evans, N. J., II, Levreault, R. M., Beckwith, S., & Skrutskie, M. 1987, ApJ, 320, 364 
Evans, N. J., II, Mundy, L. G., Kutner, M. L., & DePoy, D. L. 1989, ApJ, 346, 212 
Forster, J. R., Welch, W. J., & Wright, M. C. H. 1977, ApJ, 215, L121 

Fukui, Y., Iwata, T., Mizuno, A., Bally, J., & Lane, A. P. 1993, in Protostars and Planets III, ed. E. H. 
Levy & J. Lunine (Tucson: Univ. of Arizona Press), in press 

Genzel, R., & Downes, D. 1977, A&AS, 30, 145 

1979, A&A, 72, 234 

Georgelin, Y. M., & Georgelin, Y. P. 1976, A&A, 49, 57 

Giovanardi, C, Lizano, S., Natta, A., Evans, N. J., II, & Heiles, C. 1992, ApJ, 397, 214 
Glassgold, A. E., Mamon, G. A., & Huggins, P. J. 1989, ApJ, 336, L29 
1991, ApJ, 373, 254 

Grevasse, N., Lambert, D. L., Sauval, A. J., van Dishoeck, E. F., Farmer, C. B., & Norton, R. H. 1991, 

A&A, 242, 488 

Hackwell, J. A., Grasdalen, G. L., & Gehrz, R. D. 1982, ApJ, 252, 250 



- 22 - 



Harvey, P. M., Campbell, M. F., Hoffmann, W. F., Thronson, H. A., Jr., & Gatley, I. 1979, ApJ, 229, 990 
Harvey, P. M., & Forveille, T. 1988, A&A, 197, L19 

Hayashi, M., Hasegawa, T., Omodaka, T., Hayashi, S. S., & Miyawaki, R. 1987, ApJ, 312, 327 
Ho, P. T. P., Moran, J. M., & Rodriguez, L. F. 1982, ApJ, 262, 619 

Hollcnbach, D., & McKee, C. F. 1989, ApJ, 342, 306 

Howe, J. E., JafFe, D. T., Genzel, R., Stacey, G. J., Townes, C. H., & Lugten, J. B. 1987, BAAS, 19, 1015 
Hughes, V. A. 1985, ApJ, 298, 830 
Koo, B.-C. 1989, ApJ, 337, 318 
1990, ApJ, 361, 145 

Koppenaal, K., Sargent, A. I., Nordh, L., van Duinen, R. J., & Aalders, J. W. G. 1979, A&A, 75, LI 
Kraus, J. D. 1966, Radio Astronomy (New York: McGraw Hill) 
Kutner, M. L., & Ulich, B. L. 1981, ApJ, 250, 341 
Lada, C. J. 1985, ARA&A, 23, 267 

Lada, C. J., Blitz, L., Reid, M. J., & Moran, J. M. 1981, ApJ, 243, 769 

Lada, C. J., Thronson, H. A., Smith, H. A., Schwartz, P. R., & Glaccum, W. 1984, ApJ, 286, 302 
Langer, W. D., & Penzias, A. A. 1990, ApJ, 357, 477 

Lester, D. F., Harvey, P. M., Joy, M., & Ellis, H. B., Jr. 1986, ApJ, 309, 80 

Lizano, S., Heiles, C., Rodriguez, L. F., Koo, B.-C., Shu, F. H., Hasegawa, T., Hayashi, S., & Mirabel, I. F. 
1988, ApJ, 328, 763 

Mangum, J. G. 1993, PASP, in press 

Margulis, M., & Lada, C. J. 1985, ApJ, 299, 925 

Margulis, M., & Snell, R. L. 1989, ApJ, 343, 779 

Masson, C. R., & Chernin, L. M. 1993, ApJ, in press 

Masson, C. R., Mundy L. G., & Keene, J. 1990, ApJ, 357, L25 

McGregor, P. J., Persson, S. E., & Cohen, J. G. 1984, ApJ, 286, 609 

Merrill, K. M., & Soifer, B. T. 1974, ApJ, 189, L27 

Mitchell, G. F., & Hasegawa, T. I. 1991, ApJ, 371, L33 

MitcheU, G. F., Hasegawa, T. I., & Schella, J. 1992, ApJ, 386, 604 

Mitchell, G. F., Maillard, J.-R, & Hasegawa, T. I. 1991, ApJ, 371, 342 

Natta, A. & Giovannardi, C. 1990, ApJ, 356, 646. 

Natta, A., Giovanardi, C, Palla, F., & Evans, N. J., II 1988 ApJ, 327, 817 
Neufeld, D. A., & Dalgarno, A. 1989a, ApJ, 340, 869 
1989b, ApJ, 344, 251 

Panagia, N. 1991, in The Physics of Star Formation and Early Stellar Evolution, ed. C. J. Lada & N. D. 
Kylafis (Dordrecht: Kluwer), 565 

Parker, N. D., Padman, R., & Scott, P. F. 1991, MNRAS, 252, 442. 



- 23 - 



Persson, S. E., Geballe, T. R., Simon, T., Lonsdale, C. J., & Baas, F. 1981, ApJ, 251, L85 
Phillips, J. P., et al. 1988, A&A, 190, 289 

Plambeck, R. L., Snell, R. L., & Loren, R. B. 1983, ApJ, 266, 321 

Plume, R., Jaffc, D. T., & Evans, N. J., II 1992, ApJS, 78, 505 

Raga, A. C, Canto, J., Calvct, N., Rodriguez, & Torrelles, J. M. 1993, A&A, in press 

Reynolds, S. P. 1986, ApJ, 304, 713 

Richardson, K. J., White, G. J., Avery, L. W., Lesurf, J. C. G., & Harten, R. H. 1985, ApJ, 290, 637 
Richer, J. S., Hills, R. E., & Padman, R. 1992, MNRAS, 254, 525 
Rodriguez, L. F., & Canto, J. 1983, Rev. Mexicana Astr. Ap., 8, 163 
Rubin, R. H. 1968, ApJ, 154, 391 

Ruden, S. P., Glassgold, A. E., & Shu, F. H. 1990, ApJ, 361, 546 
Russell, A. P. G., Bally, J., Padman, R., & Hills, R. E. 1992 ApJ, 387, 219. 
Simon, M., Felli, M., Cassar, L., Fischer, J., & Massi, M. 1983, ApJ, 266, 623 
Snell, R. L., & Bally, J. 1986, ApJ, 303, 683 

Snell, R. L., Loren, R. B., & Plambeck, R. L. 1980, ApJ, 239, L17 

Snell, R. L., Scoville, N. Z., Sanders, D. B., & Erickson, N. R. 1984, ApJ, 284, 176 

Stacey, G. J. 1992, private communication 

Stacey, G. J., Jaffe, D. T., Lugten, J. B., Genzel, R., & Townes, C. H. 1987, BAAS, 19, 1016 

Stabler, S. W. 1993, in Astrophysical Jets, ed. M. Livio, C. O'Dea, & D. Burgarella (Cambridge: Cambridge 
Univ. Press), in press 

Thompson, R. I. 1984, ApJ, 283, 165 

van Dishocck, E. F., Glassgold, A. E., Guclin, M., Jaffe, D. T., Neufeld, D., Ticlcns, A. G. G. M., & 

Walmsley, C. M. 1992, in lAU Symp. 150, Astrochemistry of Cosmic Phenomena, ed. P. D. Singh 
(Dordrecht: Kluwer), 285 

Walker, C. K., Narayanan, G., Biittgenbach, T. H., Carlstrom, J. E., Keene, J., & Phillips, T. G. 1993, in 
preparation 

Werner, M. W., et al. 1980, ApJ, 242, 601 

Wood, D. O. S., & Churchwell, E. 1989, ApJS, 69, 831 

Wynn-WiUiams, C. G., Bcckhn, E. E., Forster, J. R., Matthews, K., Neugebauer, G., Welch, W. J., & 
Wright, M. C. H. 1977, ApJ, 211, L89 

Wynn-Wihiams, C. G., Beckhn, E. E., & Neugebauer, G. 1972, MNRAS, 160, 1 



This preprint was prepared with the AAS lAT^jX macros v3.0. 



-24- 



FIGURE CAPTIONS 

Fig. 1. — (a) CO spectra of W3 IRS 5 with baselines. The upper four panels show the full spectra of 
the lines, and the middle and lower four panels have blown-up vertical scales to show the details of 
the wings. All the horizontal axes cover the same velocity extent. Vertical dashed lines below the 
baselines are the wing boundaries (also given in Table 4). (6) Contour maps of the integrated intensity 
for each wing of W3 IRS 5 in the CO J = 3^2 line. The lower panel shows the maps of the inner 
wings, and the upper panel shows those of the outer wings. Solid contours are for red wings, and 
dashed contours for blue wings. At the upper right corner of each panel are the integrated antenna 
temperature (J T^dV) of the lowest contour level, the highest contour level, and the spacing between 
contours in K km s~^. Maps were made with a 12" grid. Maps also include markers at the positions 
of IRS 6, IRS 5, IRS 7, and IRS 3 (•) from left to right (Wynn- Williams, Becklin, & Neugebauer 
1972), and H2O maser sources (A) from Forster, Welch, & Wright (1977). The maps of the outer 
wings show a compact source centered at or near IRS 5. The emission in the red wing appears slightly 
extended toward the southwest (cf. Mitchell et al. 1991; Claussen ct al. 1984). (c) At the central 
position of W3 IRS 5, the line ratio (Tr) of CO J = 3^2 to J = 2^1 in the top panel, CO J = 3^2 
to ^^CO J = 3^2 in the second panel, optical depth of CO J = 3^2 line in the third panel, and line 
ratio of CO J — 2-^1 to ^^CO J = 2^1 in the bottom panel. Values for line ratios have been left 
out whenever Tjj's drop below the thresholds given in Table 3. Vertical dashed lines are the wing 
boundaries. Solid exponential curves in the inner wings in the third panel are the least-squares fits of 
T32. For the outer wings, we assumed T32 <C 1. The line ratios ^^R32, **°J?3, and **°i?2 rise toward the 
line wings until they become undefined because the line in the divisor is too weak. The ratio ^^ii32 
implies low optical depth, even for the inner wings. Even assuming the maximum ratio for the beam 
filling factors, ^^i?32 becomes large enough to require r32 < 1 for inner and outer wings, contrary to 
the T32 calculated from ^^"R^. In fact, ^^R32 is illegally high [> (3/2)^] at some channels in the outer 
wings, perhaps because of a pointing error in the CO J = 2^1 data (note that the blue EHV outflow 
of W3 IRS 5 is the most compact among the sources in this paper). Nonetheless, we interpret the 
high ratios as strong evidence for optically thin emission. 

Fig. 2. — (a) The same as Fig. la for GL 490. Only GL 490 has CO J = 2— »1 wings as strong as the 
CO J = 3— >2 wings. The J = 1^0 line is also comparable in strength to the J = 2^1 line over the 
velocity interval of the inner wings (Margulis & Lada 1985). (h) The same as Fig. lb for GL 490. 
Maps were made with a 9" grid. For the maps of red wings the lowest contours are the southernmost 
ones. Markers indicate the position of an infrared source (•) from McGregor, Persson, & Cohen 
(1984). The radio continuum source is coincident with the infrared source. Our map does not cover 
the whole outflow region and does not reveal a clear morphology, (c) The same as Fig. Ic for GL 490. 
The ratio ^^i?32 is ~0.9 throughout the line core and the wings, suggesting that the gas is cold {Tg^ 
^ 50 K) and/or the wings are optically thick, even at extremely high velocities. So we assumed that 
T32 is 2 in the EHV wings. 

Fig. 3. — (a) The same as Fig. la for NGC 2071. There is a narrow {AV f» 20 km s~-^) spectral feature 

at ^55 km near the southwestern edge of our map, consistent with the feature found by Chernin 
& Masson (1992). (6) The same as Fig. 16 for NGC 2071. Maps were made with a 9" grid. Markers 
indicate IRS 4, IRS 2, IRS 3, and IRS 1 (•) from northeast to southwest (Persson et al. 1981) and 
an H2O maser source (A) from Genzel & Downes (1979). The bipolar nature of the outflow is clear 
although our map does not completely cover the blue and red peaks (cf. Bally 1982). Chernin & 
Masson (1992) have mapped the EHV gas to larger distance and found that it covers a smaller region 
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than the HV gas. (c) The same as Fig. Ic for NGC 2071. The isotopic ratio, '*°J?25 in the inner wings 
agrees well with that of SneU ct al. (1984). The interpretation of the hne ratios is uncertain for this 
source because none of the wings have a peak at the central position. Both **°i?3 and ''*°i?2 show the 
typical rise toward the outer wings. The red side of '^^"R^ and the blue side of '*°i?2 seem to decrease 
beyond \V — Vol ~ 15 km s~^, but the uncertainties are very large. The ratio ^^i?32 rises only to the 
red side of the spectrum. 

Fig. 4. (a) The same as Fig. la for W28 A2. The peak at V = 100 km s"! in the CO J = 3^2 spectrum 
is the CS J = 7^6 line in the other sideband. Note that even the CS line is as extended as our 
inner wings, which suggests that the inner wings are emitted by dense gas (Plume, Jaffe, & Evans 
1992). The baseline of the CO J = 3^2 line is quite poor because of the broad CS line. Since the 
inner red wing of CO J = 3-^2 contains a dip which may be either emission in the reference beam or 
self- absorption, we define the inner red wing using the -^^CO J = 3— >2 line. (&) The same as Fig. lb 
for W28 A2. Maps were made with a 12" grid. The map of the inner red wing is omitted because the 
spectra contain a dip. Markers indicate H2O maser sources (A) from Genzel & Downes (1977), and a 
compact H II region (x) from Wood & Churchwell (1989). The contour map shows a compact source 
with both blue and red wings peaking at the position of an infrared source and H2O maser and also 
well within the extent of the radio continuum emission, (c) The same as Fig. Ic for W28 A2. The 
line ratios in the outer wings are somewhat unreliable because of poor baselines. The ratios **°i?2 and 
'*°i?3 rise toward the wings until they become undefined, which indicates that the optical depth is 
getting smaller as the velocity increases. 

Fig. 5.— (a) The same as Fig. la for GL 2591. The peak at F = 90 km s'^ is the CS J = 7^6 line. (6) 
The same as Fig. 16 for GL 2591. Maps were made with an 18" grid. Markers indicate an infrared 
source (•) from Lada et al. (1984), an H2O maser source (A) from Wynn- Williams et al. (1977), and 
radio continuum peaks (x) from Campbell (1984). The integrated intensities of the inner wings do 
not peak within our map. The CO J = 2^1 map of Mitchell et al. (1991) shows that the peaks of 
the inner wings are indeed beyond the boundaries of our map. Both outer wings clearly have compact 
emission peaked at the same position. In this region, the H2O maser source, the infrared source, and 
the strongest radio continuum source are all displaced from each other, lying along a NE-SW line. 
The peak of the outer wing emission is displaced from all these sources, but is closest to the radio 
continuum source. The offset from all the sources is smaller than the beam size, but larger than 
the combined pointing uncertainties and positional uncertainties of the maser and infrared source. 
If observations with higher resolution confirm this offset, it would suggest that the source of radio 
continuum emission, rather than the infrared source, drives the outflow, (c) The same as Fig. Ic 
for GL 2591. The ratios *'*°i?3 and ^^i?32 have minima at the line center and rise toward the line 
wings until they become undefined, suggesting that the outer wings are optically thin. Since only one 
channel in the inner red wing has a well-defined '*''J?3, we included the outermost channel of the line 
core to make the r32 fit. 

Fig. 6.— (a) The same as Fig. la for S140. The CO J = 3^2 spectrum shows two CS J = 7^6 lines 

(at V = —105 and 90 km s"'^) because it is a sum of data taken during two periods when the source 
had two different topocentric velocities, (b) The same as Fig. lb for S140. Maps were made with a 
12" grid. Markers indicate infrared and radio sources, IRS 2, IRS 1, and IRS 3 (•) clockwise from 
the top (Hackwell, Grasdalen, & Gehrz 1982; Evans et al. 1989), and an H2O maser source (A) from 
Genzel & Downes (1979). The inner wings show a clear bipolar pattern (cf. Hayashi et al. 1987). The 
contour map of the outer wings looks clumpy because of several secondary peaks. S140 is the only 
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source in our sample in which the outer wings appear clumpy, (c) The same as Fig. Ic for S140. The 
isotopic ratio, *''°i?2, of the inner wings is consistent with SncU et al. (1984). Both isotopic ratios rise 
toward the outer wings. The ratio ^^i?32 rises toward the blue wing, but decreases toward the red 
wing, perhaps a result of observing with different beam sizes a wing whose spatial distribution does 
not peak at the center of the source; i.e., the CO J = 2^1 beam is larger than the CO J = 3— >2 
beam and contains more of the red peak. We assumed t^2 = 2 in the red wing. 

Fig. 7.— (a) The same as Fig. la for Cepheus A. The peak at F = -110 km s"^ is CS J = 7^6 line, [h) 
The same as Fig. 16 for Cepheus A. Maps were made with a 12" grid. Markers indicate a mid-infrared 
source (•) from Beichman, Becklin, & Wynn- Williams (1979), radio component 2 of Hughes (1985) 
(x), and an H2O niaser somcc (A) from Lada ct al. (1981). The contour maps clearly show the 
bipolar nature of the inner wings. The outer wings, however, show a rather complex structure. The 
outer blue wing has a local peak at the center position, while the red wing has a displaced peak. 
The peak of the blue wing coincides within the uncertainties with an H2O maser source and a radio 
contimmm source, (c) The same as Fig. Ic for Cepheus A. The ratios ''*°J?,3 and '*°i?2 have minima at 
the line center and rise toward the line wings until the denominators become less than the thresholds. 
Since ^^i?32 is also consistent with low T32, we have assumed optically thin emission in the outer wings. 

Fig. 8. — D{n,TK) (defined by eq. [6]) is plotted for reasonable values of n (10^-10^ cm~^) and Tk 
(10-200 K), as calculated from the LVG model, and from the LTE assumption. In the LVG 

calculations, we assumed a constant CO column density per unit velocity interval, N/dV = IQ^^ cm~^ 
km^^ s, typical of our derived values, but D(n,TK) is insensitive to N/dV. 

Fig. 9. — The ratio of EHV to HV driving force versus luminosity {upper panel) and the rate of production 
of Lyman-continuum photons {lower panel). We assumed that the EHV gas is the stellar wind 
([C]/[CO] = 1) and that the HV outflow is swept-up molecular gas ([C]/[CO] = 8). Filled pointers are 

from this paper, and open pointers are from references given in Tabic 6. If there arc multiple references 
for one source, corresponding pointers are connected by vertical lines. See Table 6 for details. 

Fig. 10. — ^^i?32 as a function of T32 for various T^^- It is assumed that T21 follows eq. (4). Numbers on the 
left vertical axis are for (f) = (3/2)^ (point source), and those on the right vertical axis are for </> = 1.0 
(extended source). 

Fig. 11. — The observed excitation temperature, Tex,obsi as a function of the real excitation temperature. 
Tea;, with two components in the beam. It is assumed that X = 56, t„a — 1, TuB = 10, Tg^A ~ T^xB, 
and all the filling factors are the same. The solid curve is for the combination of CO J = 2^1 and CO 
J = 1—^0, and the dashed curve is for the combination of CO J = 3— >2 and CO J = 2^1 transitions. 
The straight dotted line is Tex,obs = Tex- 



